UNCLASSIFIED 


Defense  Technical  Information  Center 
Compilation  Part  Notice 

ADP010515 

TITLE:  A System  for  the  Aerodynamic  Optimization 
of  Three-Dimensional  Configurations 

DISTRIBUTION:  Approved  for  public  release,  distribution  unlimited 


This  paper  is  part  of  the  following  report: 

TITLE:  Aerodynamic  Design  and  Optimisation  of 
Flight  Vehicles  in  a Concurrent 
Multi-Disciplinary  Environment  [la  Conception  et 
T optimisation  aerodynamiques  des  vehicules 
aeriens  dans  un  environnement  pluridisciplinaire 
et  simultane] 

To  order  the  complete  compilation  report,  use:  ADA388284 

The  component  part  is  provided  here  to  allow  users  access  to  individually  authored  sections 
of  proceedings,  annals,  symposia,  ect.  However,  the  component  should  be  considered  within 
the  context  of  the  overall  compilation  report  and  not  as  a stand-alone  technical  report. 

The  following  component  part  numbers  comprise  the  compilation  report: 


19-1 
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Lilienthalplatz  7,  D-38108  Braunschweig,  Germany 


Summary 

The  paper  presents  a system  for  the  aerodynamic  optimization 
of  three-dimensional  configurations.  This  system  is  based  on  the 
repeated  calculation  of  the  flowfield  around  three-dimensional 
geometries  by  solving  the  Euler/Navier- Stokes  equations.  The 
basic  structure  of  the  system  and  the  incoperated  modules  are 
described.  Under  the  same  conditions  the  system  must  provide 
the  same  solutions  of  classic  aerodynamic  optimization  prob- 
lems as  given  in  literature.  So  the  function  of  the  system  is 
checked  with  the  Rhombus  airfoil  and  the  Sears -Haack  body. 
The  potential  of  the  system  is  demonstrated  with  current  aerody- 
namic optimization  problems. 
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1 Introduction 

Future  military  and  civilian  aerospace  vehicles  must  fulfil  exten- 
sive and  partially  conflicting  requirements.  Technical  aspects  like 
flying  characteristics  or  stealth  capabilities  have  to  be  considered 
as  well  as  economical  aspects  like  life-cycle  costs  or  fleet  com- 
patibilities and  ecological  aspects  like  noxious  emissions  or 
noise  regulations.  In  order  to  fulfil  these  requirements  the  close 
interaction  of  all  involved  disciplines  and  the  optimal  use  of  all 
technical  potentialities  is  necessary.  Simultaneously,  design  cycle 
times  have  to  be  reduced. 

Size  and  complexity  of  this  task  lead  to  growing  importance  of 
numerical  design  and  optimization  methods.  In  general,  these 
methods  can  be  divided  based  on  their  orientation  into  multidis- 
ciplinary methods  and  disciplinary  methods.  Multidisciplinary 
methods  are  used  mainly  in  the  predesign  stage  for  a rough  as- 
sessment of  the  vehicle  characteristics  and  a preliminary  interac- 
tion of  the  involved  disciplines.  An  overview  of  multidisciplinary 
design  and  optimization  methods  is  given  in  Ref.  1 . Disciplinary 
methods  are  used  mainly  in  the  design  stage  for  an  exact  calcula- 
tion of  the  vehicle  characteristics  and  their  improvement.  An 
overview  of  disciplinary  design  and  optimization  methods  is 
given  in  Ref.  2. 

Most  aerodynamic  design  and  optimization  methods  are  still 
based  on  linear  theory.  They  cannot  or  cannot  exactly  model  non- 
linear phenomenons  like  separations  or  vortices.  Hence,  results 
represent  only  approximations  of  real  optima.  In  order  to  im- 
prove these  approximations  the  Institute  of  Design  Aerodynam- 
ics of  the  German  Aerospace  Center  (DLR)  is  working  on  new 
methods3’5  based  on  non-linear  theory,  i.e.  the  solution  of  the 
Euler/Navier- Stokes  equations. 

Section  2 deals  with  the  fundamentals  of  such  an  aerodynamic 
optimization  method.  The  basic  structure  of  the  computational 
system  and  the  incoporated  modules  are  described.  Tests  of  the 
system  make  up  section  3.  Its  function  is  checked  with  solutions 
of  classic  aerodynamic  optimization  problems:  the  Rhombus  air- 
foil and  the  Sears-Haack  body.  Section  4 presents  applications  of 
the  system.  Aerodynamic  optimizations  of  a biconic  re-entry 
capsule,  an  SCT  wing  as  well  as  an  SCT  wing/body  configura- 
tion and  their  results  are  presented. 

2 Aerodynamic  Optimization  System 

The  Aerodynamic  Optimization  System  is  based  on  the  repeated 
calculation  of  flowfields  around  three-dimensional  geometries.  A 
mathematical  algorithm  gives  changes  of  the  respective  geome- 
try according  to  an  optimization  strategy.  These  changes  are  car- 
ried out,  the  spatial  mesh  is  created  and  the  flowfield  computed. 
The  changes  are  characterized  by  the  resulting  aerodynamic  co- 
efficients and  judged  with  an  objective  function.  Finally,  the 
mathematical  algorithm  determines  the  changes  of  the  respective 
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geometry,  which  fulfil  the  demands  formulated  through  the  ob- 
jective function  best. 

2.1  Basic  Structure 

The  basic  structure  of  the  Aerodynamic  Optimization  System 
and  the  course  of  events  is  depicted  in  Fig.  1.  Input  parameters 
are  start  values  for  the  geometric  and  aerodynamic  parameters, 
which  are  included  in  the  optimization  process,  constraints  lim- 
iting their  variations,  start  values  for  the  stepsizes,  which  are 
used  to  determine  gradients,  and  the  objective  function.  These 
parameters  have  influence  in  the  whole  system.  Further  input  pa- 
rameters are  values  for  the  geometric  and  aerodynamic  parame- 
ters, which  are  not  included  in  the  optimization  process,  and  in- 
structions for  the  generation  of  the  surface  mesh,  the  generation 
of  the  spatial  mesh  and  the  calculation  of  the  flowfield.  These 
parameters  have  influence  only  in  certain  modules. 

At  the  beginning  a geometry  corresponding  to  the  geometric  pa- 
rameters and  a surface  mesh  are  created.  For  that  the  instructions 
for  the  generation  of  the  surface  mesh  prescribe  the  number  of 
points  in  axial  and  spanwise  direction,  the  location  of  points 
near  edges  and  comers,  the  distribution  of  points  near  shocks 
and  separations  etc.  Based  on  the  surface  mesh  a spatial  mesh  is 
created.  The  instructions  for  the  generation  of  the  spatial  mesh 
prescribe  the  number  of  points  normal  to  the  surface,  the  loca- 
tion of  the  farfield,  the  distribution  of  points  etc. 

Then  the  flowfield  around  the  geometry  is  computed  and  the 
aerodynamic  coefficients  are  determined.  A numerical  method  is 
used  to  solve  the  Euler/Navier- Stokes  equations  at  every  point  of 
the  spatial  mesh.  The  instructions  for  the  calculation  of  the  flow- 
field  control  this  method.  They  give  the  number  of  iteration  cy- 
cles, the  damping  of  the  method  etc. 

In  the  end  the  aerodynamic  coefficients  are  put  into  the  objective 
function  and  the  function  is  analyzed.  A decision  on  the  continu- 
ation of  the  optimization  process  is  made.  In  case  the  gradients 
of  the  objective  function  do  indicate  the  possibility  of  improve- 
ments the  process  is  continued.  The  parameters,  which  are  in- 
cluded in  the  optimization  process,  are  changed  by  a mathemati- 
cal algorithm  according  to  an  optimization  strategy.  The 
constraints  are  checked  and  if  necessary  corrections  are  carried 
out.  The  described  steps  start  again,  i.e.  a geometry  correspond- 
ing to  the  changed  geometric  parameters  is  generated  etc.  In 
case  the  gradients  of  the  objective  function  do  not  indicate  the 
possibility  of  improvements,  even  after  changes  of  the  stepsizes, 
which  are  used  to  determine  the  gradients,  the  process  is 
stopped.  The  geometric  parameters  and  the  aerodynamic  coeffi- 
cients are  stored. 

So  the  aerodynamic  optimization  system  is  a repeated  cycle  of 
four  steps:  geometry  generation,  mesh  generation,  flowfield  cal- 
culation and  optimization.  These  steps  are  carried  out  within  the 
system  by  four  independent  modules.  Their  communication  is 
completely  based  on  data  files.  Here  advantages  regarding  easy 
maintenance  and  development  predominate  disadvantages  re- 
garding higher  computational  effort.  The  four  modules  are  de- 
scribed in  the  following  sections. 

2.2  Geometry  Generator 

The  module  for  the  geometry  generation  creates,  as  mentioned 
in  section  2.1,  a geometry  corresponding  to  the  geometric  pa- 
rameters and  a surface  mesh.  The  number  of  parameters,  which 
are  included  in  the  optimization  process,  effects  directly  the 


number  of  optimization  cycles  and  so  the  computational  effort. 
Hence,  the  module  for  the  geometry  generation  should  be  able  to 
define  a huge  variety  of  configurations  and  extensive  variations 
of  these  configurations  with  a minimum  of  geometric  parame- 
ters. The  available  CAD  systems  as  well  as  MegaCads , the  mod- 
ule for  the  mesh  generation,  cannot  guarantee  this.  Because  of 
this  a new  code,  MegaGeo,  has  been  developed  and  integrated. 

MegaGeo  composes  three-dimensional  configurations  out  of 
two-dimensional  curves  describing  elements  of  the  configura- 
tions. Figure  2 illustrates  two-dimensional  curves  describing  an 
airfoil.  Points  at  the  leading  and  trailing  edge  divide  the  airfoil 
into  upper  and  lower  side.  Further  points  divide  upper  and  lower 
side  into  segments.  The  number  and  position  of  these  points  de- 
pend on  the  complexity  of  the  airfoil  and  the  expected  flexibility 
of  the  parameterization.  The  order  of  the  segments  determines 
start  and  end  points.  A set  of  curve  functions  introduced  by  H. 
Sobieczky6  is  used  to  define  the  curve  in  the  segments.  In  princi- 
ple all  functions  are  controlled  by  four  parameters,  which  can  be 
associated  with  the  slope  and  curvature  at  the  start  and  end 
points.  But  some  functions  have  given  parameters  because  of 
their  mathematical  background  or  constraints.  The  key  number 
determing  upper  or  lower  side,  the  coordinates  of  the  start  point, 
the  key  number  determing  the  curve  function  and  the  parameters 
of  the  function  are  the  necessary  information  for  a segment. 

MegaGeo  distinguishes  two  major  elements  of  three-dimensional 
configurations:  Wing  and  body.  The  elements  forming  a wing  are 
shown  in  Fig.  3.  A certain  number  of  airfoils,  the  leading  edge, 
the  trailing  edge  and  the  twist  axis,  all  defined  by  curve  func- 
tions, together  with  information  about  the  transition  between  the 
airfoils  and  their  twist  build  up  the  basic  structure.  All  airfoils  are 
given  in  a normalized  way  and  the  leading  and  trailing  edge  are 
used  to  determine  their  magnitude.  The  elements  forming  a body 
are  shown  in  Fig.  4.  A certain  number  of  body  sections,  the  water 
line  and  the  body  axis,  all  defined  by  curve  functions,  together 
with  information  about  the  transition  between  the  sections  build 
up  the  basic  structure.  All  body  sections  are  given  in  a normal- 
ized way  and  the  water  line  is  used  to  determine  their  magnitude. 

The  surface  of  a wing,  a body  or  a wing/body  configuration  de- 
velops with  the  surface  mesh.  For  that  a distribution  of  planes  x = 
const,  is  given.  In  these  planes  the  intersection  between  wing  and 
body  results  from  a calculation.  Also  in  these  planes  the  number 
of  mesh  points  as  well  as  their  distribution  on  wing  and  body,  de- 
fined by  curve  functions  too,  is  given.  So  the  complete  geometry 
and  surface  mesh  generation  is  based  on  parameters  controlling 
curve  functions  and  can  be  carried  out  as  an  automatic  genera- 
tion. The  instructions  for  the  generation  of  the  geometry  and  the 
surface  mesh  are  more  or  less  equivalent  to  these  parameters. 

2.3  Mesh  Generator 

The  module  for  the  mesh  generation  creates,  as  mentioned  in 
section  2.1,  a spatial  mesh  based  on  the  surface  mesh  created  by 
the  module  for  the  geometry  generation.  It  has  to  guarantee  an 
automatic  generation  of  high  quality  within  the  optimization  pro- 
cess. A huge  variety  of  configurations  and  extensive  variations  of 
these  configurations  should  not  call  this  into  question.  Further- 
more, it  has  to  guarantee  an  automatic  generation  with  minimum 
computational  effort.  The  DLR  code  MegaCads1  fulfils  this  and 
has  been  integrated. 

MegaCads  enables  the  definition  of  structured  multi-block 
meshes  around  complex  three-dimensional  configurations  using 
a parametric  approach.  Basic  techniques  for  the  generation  of  ge- 
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ometries  and  different  algorithms  for  the  determination  of  sur- 
face/surface intersections  are  integrated.  Efficient  algebraic 
methods  for  two  and  three-dimensional  mesh  generation  are 
available.  In  order  to  smooth  the  algebraic  meshes  elliptic  algo- 
rithms can  be  used  allowing  different  controls  for  the  point  dis- 
tribution at  boundaries.  Support  features  for  multi -grid  and  chi- 
mera techniques  are  developed.  The  complete  interactive  mesh 
generation  is  stored  in  a protocol  file  using  a simple  script  lan- 
guage and  can  be  repeated  for  similar  configurations  or  within 
an  optimization  process  as  an  automatic  generation.  The  instruc- 
tions for  the  generation  of  the  spatial  mesh  are  more  or  less 
equivalent  to  this  protocol  file. 

2.4  Flow  Solver 

The  module  for  the  flowfield  calculation  computes,  as  men- 
tioned in  section  2. 1 , the  flowfield  around  the  geometry  and  de- 
termines the  aerodynamic  coefficients.  The  flowfield  calculation 
solving  the  Euler/Navier- Stokes  equations  has  by  far  the  greatest 
share  of  the  computational  effort.  Hence,  the  module  should  be 
able  to  carry  out  one  calculation  with  minimum  computational 
effort  as  well  as  high  robustness  and  accuracy.  The  DLR  code 
CEVCATS  and  its  successor  F LOWer 8 fulfil  this  and  have  been 
integrated. 

CEVCATS  and  FLOWer  solve  the  Euler/Navier- Stokes  equations 
using  a second-order  accurate  finite- volume  method.  The  spatial 
discretization  is  done  by  means  of  a cell- vertex  scheme  and  cen- 
tral differences.  In  order  to  damp  numerical  oscillations  first  and 
second-order  dissipative  terms  are  added  to  the  governing  equa- 
tions according  to  A.  Jameson9.  The  integration  in  time  of  the 
resulting  system  of  ordinary  differential  equations  is  carried  out 
by  means  of  an  explicit  five-stage  Runge-Kutta  scheme.  To  ac- 
celerate convergence  local  time  stepping,  residual  smoothing 
and  multi-grid  algorithms  are  used.  The  instructions  for  the  cal- 
culation of  the  flowfield  are  more  or  less  equivalent  to  the  input 
file  of  CEVCTAS  or  FLOWer . 

2.5  Optimization  Algorithm 

The  module  for  the  optimization  analyses,  as  mentioned  in  sec- 
tion 2.1,  the  objective  function,  makes  a decision  on  the  continu- 
ation of  the  optimization  process,  changes  the  parameters  and 
checks  the  constraints.  It  has  to  guarantee  a certain  determina- 
tion of  global  extreme.  Complex  structures  of  the  parameter 
space  or  numerical  fluctuations  of  the  objective  function  should 
not  call  this  into  question.  Furthermore,  it  has  to  guarantee  a de- 
termination with  minimum  computational  effort.  The  algorithm 
EXTREM  developed  by  H.G.  Jakob10  fulfils  this  and  has  been 
integrated. 

EXTREM  determines  the  changes  of  the  geometry,  which  fulfil 
the  demands  formulated  through  the  objective  function  best,  ap- 
plying finite  differences  and  parabolic  interpolation.  Figure  5 
sketches  the  principle  of  the  optimization  strategy  and  the  func- 
tion of  the  mathematical  algorithm.  The  start  point  of  the  optimi- 
zation process,  the  parameter  set  Sq,  and  the  primary  search  di- 
rection AS  are  given  as  input  data.  Search  steps  Sj  = S0  - AS  and 
S2  = S0  + AS  lead  to  the  parameter  sets  Sj  and  S2.  Values  of  the 
objective  function  at  the  parameter  sets  Sq,  Sj  and  S2  are  used  for 
the  assessment  of  the  extreme  in  the  primary  search  direction  by 
means  of  parabolic  extrapolation.  The  parameter  set  S3  results. 
At  a secondary  search  direction  perpendicular  to  the  primary 
direction  is  defined  with  an  orthogonalization  procedure.  Further 
search  steps  lead  to  the  parameter  sets  S4  and  S5.  The  assessment 


of  the  extreme  in  secondary  search  direction  by  means  of  para- 
bolic extrapolation  leads  to  the  parameter  set  S6.  For  a parameter 
space  created  by  n parameters  (n-1)  secondary  search  directions 
are  defined.  Hence,  for  a parameter  space  created  by  2 parame- 
ters a so  called  step  is  completed  with  the  parameter  set  S6.  After 
such  a step  a new  primary  search  direction  passing  through  the 
last  parameter  set  and  the  last  parameter  set  of  the  previous  step 
(or  the  start  point)  is  defined.  The  described  procedure  starts 
again,  i.e.  search  steps  lead  to  new  parameter  sets  etc.  It  contin- 
ues until  the  gradients  of  the  objective  function  do  not  indicate 
the  possibility  of  improvements,  even  after  changes  of  the  step- 
sizes. 

EXTREM  considers  constraints  of  the  parameters,  which  limit 
the  variations,  depending  on  the  occurring  violation.  In  case  the 
violation  of  a constraint  is  caused  by  a search  step,  the  step  is  ne- 
glected and  the  new  parameter  set  results  from  a step  with  half  of 
the  previous  size  in  the  opposite  direction.  In  case  the  violation 
of  a constraint  is  caused  by  a parabolic  extrapolation,  the  step  is 
also  neglected  and  the  new  parameter  set  results  from  a step  with 
a quarter  of  the  calculated  size  in  the  search  direction.  Of  course 
EXTREM  checks  the  constraints  with  the  new  parameter  sets  as 
well  and  carries  out  again  necessary  corrections  in  the  same  way. 

3 Test  of  the  Optimization  System 

Several  solutions  of  classic  aerodynamic  optimization  problems 
are  documented  in  literature,  solutions  based  on  linear  theory, 
extended  linear  theory  or  exact  theory.  Under  the  same  condi- 
tions the  aerodynamic  optimization  system  must  provide  these 
solutions  as  well.  Even  under  different  conditions  regarding  the 
flowfield  calculation  the  system  must  provide  similar  solutions. 
This  has  been  checked  with  two  classic  aerodynamic  optimiza- 
tion problems. 

3.1  Rhombus  Airfoil 

The  first  problem  is  the  symmetric  airfoil  of  given  length  and 
thickness,  which  has  minimum  drag  in  supersonic  flow.  Linear 
theory11  results  in  an  airfoil  with  straight  flanks  before  and  be- 
hind the  maximum  thickness,  a so-called  rhombus  airfoil.  The 
maximum  thickness  of  this  airfoil  lies  at  x/c  = 0.5 , i.e.  exactly  at 
half  chord  length.  It  depends  not  on  the  Mach  number.  The  drag 
coefficient  amounts  to  CD  = 0.00208  for  Mach  number  M ^ = 2.0 
and  Q>  = 0.00321  for  Mach  number  M ^ = 1.5  assuming  a rela- 
tive thickness  of  b/c  = 0.03 . Extended  linear  and  exact  theory12 
result  also  in  an  airfoil  with  straight  flanks  before  and  behind  the 
maximum  thickness.  But  the  maximum  thickness  lies  not  exactly 
at  half  chord  length.  It  depends  on  the  Mach  number  and  lies  at 
x/c  = 0.5191  for  Mach  number  M «.  = 2.0  and  at  x/c  = 0.5192  for 
Mach  number  Mx  = 1.5,  i.e.  behind  half  chord  length.  The  maxi- 
mum thickness  moves  closer  to  the  trailing  edge  with  increasing 
Mach  number.  The  drag  coefficient  amounts  to  CD  = 0.00205  for 
Mach  number  Mx  = 2.0  and  CD  - 0.  00317  for  Mach  number  M ^ 
= 1.5  assuming  a relative  thickness  of  b/c  = 0.03. 

Several  optimizations  of  a symmetric  airfoil  with  respect  to  min- 
imum drag  in  supersonic  flow  and  based  on  solutions  of  the  Euler 
equations  have  been  carried  out.  They  included  five  geometric 
parameters  sketched  in  Fig.  6:  the  position  of  the  maximum 
thickness  in  x-direction,  the  slope  of  the  flank  before  the  maxi- 
mum thickness  at  the  leading  edge  and  at  the  maximum  thickness 
as  well  as  the  slope  of  the  flank  behind  the  maximum  thickness  at 
the  maximum  thickness  and  at  the  trailing  edge.  Figures  7 and  8 
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depict  initial  and  resulting  airfoils  for  Mach  number  M x = 2.0 
using  different  start  values  for  the  parameters  and  stepsizes.  All 
optimizations  lead  to  the  same  airfoil.  Hence,  the  resulting  air- 
foil is  not  a coincidence.  An  effect  of  the  start  values  could  not 
be  found.  The  optimizations  lead  to  an  airfoil  agreeing  very 
good  with  the  solution  based  on  extended  linear  and  exact  the- 
ory. The  maximum  thickness  lies  at  x/c  = 0.5157  compared  to 
x/c  - 0.5192  based  on  extended  linear  and  exact  theory.  Differ- 
ence between  both  values  is  about  1%.  The  drag  coefficient 
amounts  to  CD  = 0.00209  compared  to  CD  = 0.00205  based  on 
extended  linear  and  exact  theory.  Difference  is  about  2%. 

Figure  9 depicts  initial  and  resulting  airfoils  for  Mach  number 
Mk  = 1.5  Again,  the  optimization  leads  to  an  airfoil  agreeing 
very  good  with  the  solution  based  on  extended  linear  and  exact 
theory.  The  maximum  thickness  lies  at  x/c  = 0.5162  compared 
to  x/c  = 0.5192  based  on  extended  linear  and  exact  theory.  Dif- 
ference between  both  values  is  about  0.5%.  The  drag  coefficient 
amounts  to  CD  = 0.00322  compared  to  CD  = 0.00317  based  on 
extended  linear  and  exact  theory.  Difference  is  2%. 

Besides  the  effect  of  start  values  for  the  parameters  and  stepsizes 
the  effect  of  accuracy  as  well  as  mesh  density  has  been  investi- 
gated. As  mentioned  in  section  2.1,  a numerical  method  is  used 
to  solve  the  Euler/Navier- Stokes  equations  at  every  point  of  the 
spatial  mesh  and  to  determine  the  aerodynamic  coefficients.  The 
accuracy,  where  this  method  stops,  can  be  varied.  Resulting  air- 
foils for  Mach  number  M ^ = 2.0  determing  the  aerodynamic  co- 
efficients with  different  accuracies  are  shown  in  Fig.  10.  Up  to 
10‘4  an  effect  of  the  accuracy  cannot  be  found.  Figure  1 1 plots 
the  computational  effort  on  a NEC  SX4  for  the  optimizations 
with  different  accuracies.  Here  an  effect  can  be  found.  In  accor- 
dance to  expectations  lower  accuracies  lead  to  lower  computa- 
tional effort.  A spatial  mesh  is  necessary  solve  the  Euler  equa- 
tions and  to  determine  the  aerodynamic  coefficients.  The 
number  of  mesh  points  can  be  varied.  Resulting  airfoils  for 
Mach  number  = 2.0  using  spatial  meshes  of  different  densi- 
ties are  shown  in  Fig.  12.  Up  to  1275  points  only  a small  effect 
of  the  density  can  be  found.  Figure  13  plots  the  computational 
effort  for  the  optimizations  with  different  mesh  densities.  Again, 
in  accordance  to  expectations  lower  densities  lead  to  lower  com- 
putational effort. 

3.2  Sears-Haack  Body 

The  second  classic  aerodynamic  optimization  problem  is  the 
body  of  revolution  of  given  length  and  volume,  which  has  mini- 
mum drag  in  supersonic  flow.  Linear  theory13  results  in  body 
with  the  maximum  thickness  at  x/lB  = 0.5 , i.e.  exactly  at  half 
body  length.  The  position  depends  not  on  the  Mach  number.  The 
drag  coefficient  amounts  to  cd  ~ o .1104  for  Mach  number  M « 
= 2.0  as  well  as  for  Mach  number  M ^ = 1.5  assuming  a volume 
of  Vb/13ref  = 0.0046.  The  coefficient  depends  also  not  on  the 
Mach  number.  Extended  linear  theory14  results  in  a body  with 
the  maximum  thickness  not  at  half  body  length.  The  position  de- 
pends on  the  Mach  number  and  lies  at  x/lB  = 0.5107  for  Mach 
number  = 2.0  and  at  x/lB  = 0.5110  for  Mach  number  — 
1.5 , i.e.  behind  half  body  length.  It  moves  closer  to  the  tail  with 
increasing  Mach  number.  A drag  coefficient  cannot  be  calcu- 
lated. 

Several  optimizations  of  a body  of  revolution  with  respect  to 
minimum  drag  in  supersonic  flow  and  based  on  solutions  of  the 
Euler  equations  have  been  carried  out.  They  included  six  geo- 
metric parameters  sketched  in  Fig.  14:  the  position  of  the  maxi- 


mum thickness  in  x-  and  z-direction,  the  slope  of  the  contour  at 
the  nose,  before  the  maximum  thickness,  behind  the  maximum 
thickness  and  at  the  tail.  Some  of  these  parameters  are  given  be- 
cause of  constraints.  The  contour  has  to  be  continual  at  the  maxi- 
mum thickness.  Hence,  the  slope  of  the  contour  before  and  be- 
hind the  maximum  thickness  must  be  equal,  i.e.  zero.  The 
volume  has  to  be  constant.  Therefore,  the  slope  at  the  tail  must  be 
chosen  depending  on  the  position  of  the  maximum  thickness  and 
the  slope  at  the  nose.  Figures  15  and  16  depict  initial  and  result- 
ing bodies  for  Mach  number  M « = 2.0  using  different  start  val- 
ues for  the  parameters  and  stepsizes.  All  optimizations  lead  to 
the  same  body.  Hence,  the  resulting  body  is  not  a coincidence. 
An  effect  of  the  start  values  cannot  be  found.  The  optimizations 
lead  to  a body  agreeing  very  good  with  the  solution  based  on  ex- 
tended linear  theory.  The  maximum  thickness  lies  at  x/lB  = 
0.5028  compared  to  x/lB  = 0.5107  based  on  extended  linear  the- 
ory. Difference  between  both  values  is  about  1.5%.  The  drag  co- 
efficient amounts  to  Cp  = 0.0889. 

Figure  17  depicts  initial  and  resulting  bodies  for  Mach  number 
M = 1.5.  Again,  the  optimization  leads  to  a body  agreeing  very 
good  with  the  solution  based  on  extended  linear  theory.  The  max- 
imum thickness  lies  at  x/lB  = 0.4942  compared  to  x/lB  = 0.5110 
based  on  extended  linear  theory.  Difference  between  both  values 
is  about  3%.  The  drag  coefficient  amounts  to  Cp  = 0.0949. 

Resulting  bodies  for  Mach  number  Mx  = 2.0  determing  the  aero- 
dynamic coefficients  with  different  accuracies  are  shown  in  Fig. 
18.  Up  to  1 0'4  an  effect  of  the  accuracy  cannot  be  found.  Figure 
19  plots  the  computational  effort  on  a NEC  SX4  for  the  optimi- 
zations with  different  accuracies.  Here  an  effect  can  be  found.  In 
accordance  to  expectations  lower  accuracies  lead  to  lower  com- 
putational effort.  Resulting  bodies  for  Mach  number  Mx  = 2.0 
using  spatial  meshes  of  different  densities  are  shown  in  Fig.  20. 
Up  to  14985  points  only  a small  effect  of  the  density  can  be 
found.  Figure  21  plots  the  computational  effort  for  the  optimiza- 
tions with  different  mesh  densities.  Again,  in  accordance  to  ex- 
pectations lower  densities  lead  to  lower  computational  effort. 

So  the  aerodynamic  optimization  system  provides  the  same  solu- 
tions for  classic  aerodynamic  optimization  problems  as  given  in 
literature.  Input  parameters,  accuracy  of  the  flowfield  calculation 
and  density  of  the  spatial  mesh  have  no  effect  on  the  solution. 
But  they  have  a significant  effect  on  the  computational  effort. 

4 Application  of  the  Optimization  System 

Different  kinds  of  aerodynamic  optimization  problems  are  en- 
countered during  the  design  of  aerospace  vehicles.  With  several 
examples  the  potential  of  the  Aerodynamic  Optimization  System 
has  been  proven.  Three  examples  are  documented  in  the  follow- 
ing sections. 

4.1  Re-entry  Capsule 

The  aerodynamic  characteristics  of  re-entry  vehicles  are  very  im- 
portant for  their  mission.  They  decide  the  time  frame  of  the  re- 
entry, the  place  for  the  landing,  the  strain  on  the  crew  as  well  as 
the  necessity  for  reactive  or  aerodynamic  control  systems. 
Hence,  an  improvement  of  the  characteristics  increases  the  possi- 
bilities for  the  mission  and  occasionally  the  financial  require- 
ments. 

Several  optimizations  of  a biconic  re-entry  capsule  with  respect 
to  various  objective  functions  and  based  on  solutions  of  the  Euler 
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equations  have  been  carried  out.  The  included  seven  geometric 
parameters  are  sketched  in  Fig.  22:  the  radius  of  the  nose,  the 
angles  of  the  first  and  second  cone,  the  length  of  the  first  and 
second  cone,  the  inclination  angle  of  the  first  cone  and  the  diam- 
eter of  the  base.  Some  of  these  parameters  are  given  because  of 
constraints,  which  have  been  taken  into  account  for  a realistic 
background.  The  radius  of  the  nose  has  to  be  larger  than 
rCf/^REF  = in  order  to  guarantee  reasonable  heat  loads.  The 
diameter  of  the  base  has  to  be  dCB/lREF  = 0.35  for  good  struc- 
tural contact  with  launch  systems.  Finally,  the  volume  of  the 
capsule  has  to  be  V(/l3REF  = 0.03  in  oder  to  accommodate  crew, 
systems  etc.  Hence,  the  length  of  the  second  cone  must  be  cho- 
sen depending  on  the  angles  of  the  first  and  second  cone  and  the 
length  of  the  first  cone. 

Figure  23  depicts  resulting  capsules  for  Mach  number  = 6.0 
with  and  without  taking  the  stability  into  account.  The  objective 
function  is  identical.  Taking  the  stability  not  into  account  leads 
to  a long  and  slender  capsule.  A difference  between  the  first  and 
second  cone  angle  cannot  be  found.  Taking  the  stability  into  ac- 
count by  moving  the  center  of  gravity  leads  to  a shorter  capsule. 
A difference  between  the  first  an  second  cone  angle  can  be 
found.  Lift-to-drag  ratios  depending  on  the  angle  of  attack  are 
plotted  in  Fig.  24.  Taking  the  stability  not  into  account  causes 
significantly  higher  lift-to-drag  ratios. 

Figure  25  depicts  resulting  capsules  for  Mach  number  M ^ = 6.0 
using  different  objective  functions.  The  stability  is  taken  into  ac- 
count. High  weighting  on  the  lift-to-drag  ratio  leads  to  a long 
and  slender  capsule  again.  High  weighting  on  the  volume  effi- 
ciency leads  to  a shorter  capsule.  As  expected  high  weighting  on 
the  lift-to-drag  ratio  causes  higher  lift-to-drag  ratios,  see  Fig.  26. 
In  this  case  bending  of  the  capsule  by  an  inclination  angle  of  the 
first  cone  increases  the  ratios  even  more.  In  the  other  case,  high 
weighting  on  the  volume  efficiency,  bending  provides  no  benefit. 

Figure  27  depicts  resulting  capsules  for  Mach  number  -6.0 
and  Mach  number  M « = 1.5.  The  objective  function  is  identical, 
the  stability  is  taken  into  account.  For  Mach  number  -1.5 
the  resulting  capsule  has  a greater  radius  of  the  nose.  Its  lift-to- 
drag  ratios  are  lower,  see  Fig.  28. 

4.2  SCT  Wing 

The  aerodynamic  characteristics  of  the  wing  are  crucial  for  an 
aircraft.  The  wing  produces  by  far  the  greatest  share  of  the  drag. 
Hence,  an  improvement  of  the  characteristics  decreases  the  fuel 
consumption  and  by  this  way  increases  the  economic  efficiency 
of  the  aircraft. 

Several  optimizations  of  the  wing  for  a Supersonic  Commercial 
Transport  with  respect  to  minimum  drag  at  cruise  conditions  and 
based  on  solutions  of  the  Euler  equations  have  been  carried  out. 
They  included  31  geometric  parameters.  The  parameters  con- 
cerning the  planform  of  the  wing  are  sketched  in  Fig.  29:  the  po- 
sition of  the  points  at  the  wing  tip  in  x-direction  and  in  x-  as  well 
as  y-direction,  the  position  of  the  points  at  the  leading-edge  kink 
and  the  trailing-edge  kink  in  x-  as  well  as  y-direction  and  the  po- 
sition of  the  point  at  the  trailing  edge  root  in  x-direction.  All 
points  are  connected  with  straight  lines.  The  parameters  con- 
cerning the  airfoils  at  the  root,  at  half  of  the  half-span  and  at  the 
tip  are  similar  to  those  in  Fig.  6:  the  position  of  the  maximum 
thickness  in  x-direction,  the  slope  of  the  flank  before  the  maxi- 
mum thickness  at  the  leading  edge  and  at  the  maximum  thick- 
ness, the  slope  of  the  flank  behind  the  maximum  thickness  at  the 
maximum  thickness  and  at  the  trailing  edge,  the  shift  of  the  air- 


foils in  z-direction,  the  chamber  and  the  local  angle  of  attack. 
Some  of  these  parameters  are  given  because  of  constraints, 
which  have  been  taken  into  account  for  a realistic  background. 
The  area  of  the  wing  has  to  be  Assw/l2REF  = 4.5.  Hence,  the  posi- 
tion of  the  trailing  edge  must  be  chosen  depending  on  the  span  of 
the  wing  as  well  as  the  position  of  the  leading-edge  kink  and  the 
trailing-edge  kink.  The  contour  has  to  be  continual  at  the  maxi- 
mum thickness  of  the  airfoils.  As  a result  the  slop  of  the  flanks 
before  and  behind  the  maximum  thickness  must  be  equal,  i.e. 
zero.  The  volume  of  the  wing  has  to  be  Vssv/l3 REF  = 0.0167  and 
the  maximum  thickness  of  the  airfoils  has  to  be  §/c  = 0.025. 
Therefore,  the  slopes  at  the  trailing  edge  must  be  chosen  depend- 
ing on  the  positions  of  the  maximum  thickness  and  the  slopes  at 
the  leading  edge.  Furthermore,  the  lift  coefficient  of  the  wing  has 
to  be  CE  = 0.12  at  Mach  number  = 2.0.  This  is  ensured  by 
the  objective  function. 

Figure  30  depicts  the  initial  and  the  resulting  wing  of  one  optimi- 
zation. The  resulting  wing  has  a greater  span  than  the  initial 
wing.  In  order  to  keep  the  area  its  length  is  smaller.  The  resulting 
wing  has  a nearly  straight  leading  edge.  The  kink  has  vanished. 
Figure  3 1 depicts  the  initial  and  the  resulting  airfoils  in  a normal- 
ized way.  Chamber,  local  angle  of  attack  etc.  are  very  small  and 
have  been  neglected  here.  The  resulting  airfoils  are  similar.  They 
have  a position  of  the  maximum  thickness  slightly  before  half 
chord  length.  The  lift-to-drag  ratio  of  the  resulting  wing  amounts 
at  = 2.0  to  L/D  = 16. 

4.3  SCT  Wing/Body  Configuration 

Several  optimizations  of  an  SCT  wing/body  configuration  with 
respect  to  minimum  drag  at  cruise  conditions  and  based  on  solu- 
tions of  the  Euler  equations  have  also  been  carried  out.  The  in- 
cluded nine  geometric  parameters  are  sketched  in  Fig.  32:  the 
length  and  diameter  of  the  cylindrical  body  part,  the  length  and 
diameter  of  the  nose  and  tail  cone,  the  position  of  the  wing  in  x- 
and  z-direction  and  the  inclination  angle  of  the  body.  Nose  and 
tail  cone  have  been  shaped  like  Sears-Haack  bodies,  the  wing  has 
been  taken  from  Ref.  15  and  kept  constant.  Some  of  the  parame- 
ters are  given  because  of  constraints,  which  have  been  taken  into 
account  for  a realistic  background.  The  body  has  to  accommo- 
date 250  passengers  at  6 seats  per  row.  Hence,  according  to 
guidelines  the  body  must  have  a diameter  of  dSSB/lREF  - 0.35  for 
6 seats  per  row  and  based  on  this  a length  of  Iss^REF  ~ $*71  for 
250  passengers.  The  contour  has  to  be  continual  at  the  transition 
from  the  nose  cone  to  the  cylindrical  body  part  and  from  the  cy- 
lindrical body  part  to  the  tail  cone.  Therefore,  the  nose  and  tail 
cone  must  have  a diameter  of  ds$(VlREF  = 0.35.  Similar  to  the  op- 
timization of  the  SCT  wing  the  lift  coefficient  of  the  configura- 
tion has  to  be  CE  = 0.1  at  Mach  number  M ^ = 2.0  and  angle  of 
attack  a = 3°.  This  is  ensured  by  the  objective  function.  An  as- 
sessment of  viscous  drag  by  the  flat  plate  analogy  is  included  in 
the  objective  function. 

Fig.  33  depicts  the  initial  and  the  resulting  configurationof  an  op- 
timization. The  resulting  configuration  has  a significantly  greater 
length  of  the  nose  and  tail  cone  than  the  initial  configuration.  It 
has  a position  of  the  wing  closer  to  the  nose  and  a greater  inclina- 
tion angle  of  the  body.  Lift-to-drag  ratios  depending  on  the  angle 
of  attack  are  plotted  in  Fig.  34.  The  resulting  configuration  has 
significantly  higher  lift-to-drag  ratios  than  the  initial  configura- 
tion, not  only  in  the  design  point,  i.e.  at  angle  of  attack  a = 3°, 
but  also  at  other  angles  of  attack. 
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5 Conclusions 

The  described  Aerodynamic  Optimization  System  has  proven  to 
be  a useful  tool  for  the  aerodynamic  optimization  of  three-di- 
mensional configurations.  The  system  is  based  on  the  repeated 
calculation  of  the  fiowfield  around  three-dimensional  geometries 
by  solving  the  Euler/Navier- Stokes  equations.  Hence,  non-linear 
phenomenons  influence  the  optima.  Under  the  same  conditions 
the  system  has  provided  the  same  solutions  of  classic  aerody- 
namic optimization  problems  as  given  in  literature.  Input  param- 
eters, accuracy  of  the  flowfield  calculation  and  density  of  the 
spatial  mesh  had  in  a reasonable  range  no  effect  on  the  solutions. 
But  they  had  a significant  effect  on  the  computational  effort.  The 
Aerodynamic  Optimization  System  has  shown  its  potential  for 
current  aerodynamic  optimization  problems.  Improvements  of  a 
re-entry  capsule,  an  SCT  wing  and  an  SCT  wing/body  configu- 
ration before  a realistic  background  have  been  achieved. 
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Figure  1 : Basic  structure  of  the  Aerodynamic  Optimization  System. 
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O resulting  airfoil  exact  theory 

Figure  7:  Initial  and  resulting  airfoils  of  optimizations  at 
Mach  number  M ^ = 2.0  using  different  start 
values  for  the  parameters. 
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resulting  airfoil  accuracy  10'4  resulting  airfoil  accuracy  10’8 

resulting  airfoil  accuracy  10’  resulting  airfoil  accuracy  lO*10 

Figure  10:  Resulting  airfoils  of  optimizations  at  Mach 
number  Mx  = 2.0  determing  the  aerody- 
namic coefficients  with  different  accuracies. 
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Figure  8:  Resulting  airfoils  of  optimizations  at  Mach 
number  M -2.0  using  different  start  values 
for  the  stepsizes. 
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Figure  9:  Initial  and  resulting  airfoils  of  an  optimizations  at 
Mach  number  = 7.5. 


Figure  11:  Computational  effort  for  the  optimizations  at 
Mach  number  M ^ = 2.0  determing  the  aero- 
dynamic coefficients  with  different  accuracies, 
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resulting  airfoil  19705  points 


Figure  1 2:  Resulting  airfoils  of  optimizations  at  Mach 
number  MK  = 2.0  using  spatial  meshes  of 
different  densities. 
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Figure  13:  Computational  effort  for  the  optimizations  at 
Mach  number  M = 2.0  using  spatial  meshes 
of  different  densities. 


Figure  14:  Body  of  revolution  with  6 geometric  parameters, 
volume  Vft/l3 REF  = 0.0046. 


coordinate  x/lB  — ►- 

initial  body  opt.  1 resulting  body  opt.  1 

initial  body  opt.  2 resulting  body  opt.  2 

initial  body  opt.  3 resulting  body  opt.  3 

O resulting  body  extended  linear  theory 

Figure  15:  Initial  and  resulting  bodies  of  optimizations  at 
Mach  number  Mx  = 2.0  using  different  start 
values  for  the  parameters. 


coordinate  x/L 


resulting  body  opt.  4 resulting  body  opt.  5 

resulting  body  opt.  6 

O resulting  body  extended  linear  theory 

Figure  16:  Resulting  bodies  of  optimizations  at  Mach 
number  M ^ = 2.0  using  different  start  values 
for  the  stepsizes. 


initial  body  opt.  8 resulting  body  opt.  8 

O resulting  body  extended  linear  theory 


Figure  17:  Initial  and  resulting  body  of  an  optimization  at 
Mach  number  = 1.5 . 


resulting  body  accuracy  10"4  resulting  body  accuracy  10*5 

resulting  body  accuracy  10’6 

Figure  18:  Resulting  bodies  of  optimizations  at  Mach 
number  = 2.0  determing  the  aerody- 
namic coefficients  with  different  accuracies. 
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Figure  19:  Computational  effort  for  the  optimizations  at 
Mach  number  M ^ = 2,0  determing  the  aero- 
dynamic coefficients  with  different  accuracies. 


Figure  22:  Biconic  re-entry  capsule  with  7 geometric 
parameters,  volume  V(/l*  ref  ~ 0.03,  base 
diameter  ^cb^ref  ~ 0.35. 


resulting  body  14985  points  resulting  body  109707  points 

resulting  body  828755  points 

Figure  20:  Resulting  bodies  of  optimizations  at  Mach 
number  M ^ = 2.0  using  spatial  meshes  of 
different  densities. 


Figure  21:  Computational  effort  for  the  optimizations  at 
Mach  number  = 2.0  using  spatial  meshes 
of  different  densities. 


resulting  capsule  G = 0.5  LTD  + 0.4  U,  with  stability 

resulting  capsule  G = 0.5  L/D  + 0.4  U,  without  stability 

Figure  23:  Resulting  biconic  re-entry  capsules  of  optimizations 
at  Mach  number  = 6.0  with  and  without  taking 
the  stability  into  account. 


resulting  capsule  G = 0.5  ITD  + 0.4  U,  with  stability 

resulting  capsule  G = 0.5  LTD  + 0.4  U,  without  stability 

Figure  24:  Lift-to-drag  ratios  of  the  biconic  re-entry  capsules 
resulting  from  optimizations  at  Mach  number 
= 6.0  with  and  without  taking  the  stability  into 
account. 
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resulting  capsule  G = 0.7  L/D  + 0.2  U 

resulting  capsule  G = 0.7  L/D  + 0.2  U,  bent 

resulting  capsule  G = 0.2  L/D  + 0.9  U 

resulting  capsule  G = 0.2  L/D  + 0.9  U,  bent 

Figure  25:  Resulting  biconic  re-entry  capsules  of  optimizations 
at  Mach  number  M ^ = 6.0  with  different  objective 
functions. 


resulting  capsule  G = 0.2  L/D  + 0.9  U,  M = 6.0 

resulting  capsule  G - 0.2  L/D  + 0.9  U,  M=1.5 

Figure  27:  Resulting  biconic  re-entry  capsules  of  optimizations 
at  Mach  number  M ^ = 6.0  and  at  Mach  number 
AC  = 1.5. 


resulting  capsule  G = 0.7  L/D  + 0.2  U 

resulting  capsule  G = 0.7  L/D  + 0.2  U,  bent 

resulting  capsule  G = 0.2  L/D  + 0.9  U 

resulting  capsule  G = 0.2  L/D  + 0.9  U,  bent 

Figure  26:  Lift-to-drag  ratios  of  the  biconic  re-entry  capsules 
resulting  from  optimizations  at  Mach  number 
M = 6.0  with  different  objective  functions. 


resulting  capsule  G = 0.2  L/D  + 0.9  U,  M = 6.0 

resulting  capsule  G = 0.2  L/D  + 0.9  U,  M = 1.5 

Figure  28:  Lift-to-drag  ratios  of  the  biconic  re-entry  capsules 
resulting  from  optimizations  at  Mach  number 
AC  = 6.0  and  at  Mach  number  AC  = i.5. 


Figure  29:  Planform  of  the  wing  for  a Supersonic  Commercial  Transport  with  8 
parameters,  volume  Vssw^ref  " 0.0167,  area  A ssw^REF  ~ 4.5 
and  lift  coefficient  CL  = 0.12. 
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Figure  30:  Planform  of  the  initial  and  resulting  SCT  wing  of  an  optimization  at 
Mach  number  - 2.0. 
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Figure  31:  Normalized  airfoils  of  the  initial  and  resulting 
SCT  wing  of  an  optimization  at  Mach  number 
= 2.0. 


Figure  33:  Initial  and  resulting  SCT  wing/body  configurations 
resulting  from  an  optimization  at  Mach  number 
M oc  = 2.0  and  angle  of  attack  a = 3°. 


Figure  32:  Wing/body  configuration  for  a Supersonic  Commer- 
cial Transport  with  9 parameters,  250  passengers  in  6 
seats  per  row,  wing  fixed  and  lift  coefficient  CL  = 0.1. 


Figure  34:  Lift-to-drag  ratios  of  the  initial  and  resulting 
SCT  wing/body  configuration  resulting  from 
an  optimization  at  Mach  number  Mx  = 2.0 
and  angle  of  attack  a = 3°. 
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DISCUSSION 


Session  III,  Paper  #19 


Prof  Knight  (Rutgers  University,  USA)  asked  the  largest  number  of  design  variables  the  author 
had  considered.  He  wanted  to  know  how  the  author’s  approach  would  “scale  up”  to  problems 
involving  large  numbers  of  design  variables. 

Mr  Orlowski  replied  that  the  largest  number  of  variables  had  been  35  optimizing  an  SCT 
wing.  He  noted  that  his  approach  using  Euler  equations  was  not  cheap  in  terms  of 
computational  effort.  He  believed  that  with  available  high  speed  computing  and 
experience  with  the  system  its  use  was  acceptable,  at  least  for  the  moment.  He  felt  that 
further  developments  in  computer  systems  available  to  industry  should  make  the  system 
acceptable  to  them  also. 

Mr  Doherty  (DERA,  UK)  asked  whether  it  could  be  confirmed  that  the  use  of  finite  difference 
gradient  generation  rendered  computation  time  independent  of  the  number  of  constraints,  unlike 
in  an  adjoint  approach  for  gradient  calculation. 

Mr  Orlowski  replied  that  the  constraints  affect  the  computational  effort  by  reducing  the 
number  of  parameters  that  are  really  included  in  the  optimization  process.  However,  the 
process  itself  was  independent  of  this  effect. 


